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Abstract 

The strange-anticharmed Pentaquark is a uudcs or uddcs five-quark baryon 
that is expected to be either a narrow resonance, or possibly even stable 
against strong and electromagnetic decay. We describe this hyperon here, its 
structure, binding energy and lifetime, resonance width, production mecha- 
nisms, production cross sections, and decay modes. We describe techniques 
to reduce backgrounds in search experiments and to optimize the conditions 
for Pentaquark observation. Possibilities for enhancing the signal over back- 
ground in Pentaquark searches are investigated by examining predictions for 
detailed momentum and angular distributions in multiparticle final states. 
General model-independent predictions are presented as well as those from 
two models: a loosely bound D~N molecule and a strongly-bound five-quark 
system. Fermilab E791 data, currently being analyzed, may have marginal 
statistics for showing definitive signals. Future experiments in the spirit of 
the recent CHARM2000 workshop, such as FNAL E781 and CERN CHEOPS 
with 10® — 10'' reconstructed charmed baryon events, should have sensitivity 
to determine whether or not the Pentaquark exists. 



1. Introduction 



Ordinary hadrons are mesons or baryons, whose quantum numbers can be de- 
scribed by quark-antiquark or three-quark configurations. Unusual hadrons 
that do not fit this picture would constitute new forms of hadronic matter 
- exotic hadrons. Such hadrons may have valence multiquark configurations 
such as qqqq and qqqqq. Exotic hadrons can have anomalous quantum num- 
bers not accessible to three-quark or quark-antiquark structures (open exotic 
states) or even usual quantum numbers (cryptoexotic states). Cryptoexotic 
hadrons can be identified only by their unusual dynamical properties (anoma- 
lously narrow decay widths, anomalous branching ratios, etc.). The discovery 
of exotic hadrons would have far-reaching consequences for quantum chromo- 
dynamics, for the concept of confinement, and for specific models of hadron 
structure (lattice, string and bag models). Detailed discussions of exotic 
hadron physics can be found in recent reviews [|l], H, ||, H, If. 

We consider here possible exotic hadronic states which contain quarks 
with four different flavors. The states have quark conflgurations Qqqq and 
qqqqQ, with one heavy quark Q(c,b) and also lighter quarks q(u,d,s). Their 
properties follow from the general hypothesis of "flavor antisymmetry" , by 
which quark systems characterized by the maximum possible antisymmetry 
of quark flavors (both quarks and antiquarks) are the most strongly bound. 
For instance,this means that the uuds system with ordinary ds flavors would 
be more bound than the uuds one with exotic values of charge Qc = +2, etc. 
Jaffe predicted in this spirit that for dibaryons with six light quarks, the 
most bound is the Hexaquark H = [uuddss] combination, for which not more 
than two quarks are in states with identical flavors. 

Thus, for the mesons and baryons with three types of light quark con- 
stituents (u,d,s), the "flavor antisymmetry" principle predicts that the states 
with open exotic charges and flavors are not so strongly bound and may have 
decay widths too large to be really observable. This may explain why the 
main candidates for light quark exotics have cryptoexotic characteristics or 
are states with exotic J^*" values, rather than having open exotic flavors such 
as mesons with Qc = ±2, or S* = ±2, baryons with 5* > 0, or Qc > 2 and so 
on. 

The situation can change for the exotic hadrons with four different quark 
constituents u, d, s, c or u, d, s, b, etc. For these states the flavor antisymme- 
try principle allows the existence of strongly bound states with open exotic 
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charges and flavors. Thus, Lipkin and Gignoux et al. |^ showed that 
5-quark "anticharmed" baryons (Pentaquarks) of the P° = [uudcs] and P~ 
= [uddcs] type, or analogous "anti-beauty" baryons, are the most bound in 
the 5-quark sector. There are also predictions 0, |1^ for the most bound 
tetraquark exotic meson, the Fs=[csud\. The present report focusses on such 
exotic states, with one heavy quark. At the CHARM2000 workshop a 
shorter version of the present work was presented [0. The properties of ex- 
otic pentaquark baryons [qqqqq, sqqqq, ssqqq, sssqq, ssssq, qqqss, qqsss] have 
been also discussed by Kaidalov, Grigoryan, Ferrer, Strakovsky, and others 
ll, K^ (see also and references therein). 



2. Pentaquark Binding Energy 

A very interesting situation can be realized, if there are exotic hadrons 
with heavy quarks, which are the bound states of known (quasi) stable 
particles, with masses which are below the threshold for strong decay. 
For example, if the Pentaquark is a bound state {ND~) with the mass 
M{P) < M{D~) + M{N), such a state would decay only via weak interac- 
tions and would be quasistable. On the other hand, a resonance {ND~)-sta.te 
with a mass above the strong threshold, would be a short-lived state that 
decays strongly as P ^ -|- D~ . 

The binding potential of a system is given by the difference between the 
Color Hyperfine (CH) interaction in the system and in the lightest color- 
singlet combination of quarks into which it can be decomposed. The wave 
function of the H may be written as: 

■^H = «1^6g + /5l^(AA) + 7l^(S-S+) + 5i^(H-p). (l) 

The lightest color singlet combination is the AA system at 2231 MeV. The CH 
contribution to the binding energy of the H is about 150 MeV in simple 
models of the CH interaction. Similarly, the and P~ wave functions can 
be written as: 

*po = a2^5g + /?2^(D7p) + 12^{S+D-) + ^2^(AS0), (2) 
= a3^5g + /33^(D7„) + 73^(E-D0) + ^S^CAD")- (3) 

Here, the lightest color singlet is the D~N system at 2907 MeV. The CH 
contribution to the mass splitting M{D-p) - M(P°) is the same as for the H 
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particle, again in simple models of the color hyperfine interaction [|], |^ . The 
ant i- Pent aquark is defined similarly. In general, whatever can be said about 
the Pentaquark holds true also for the charge-conjugate particles. 



The calculations of Ref. [|T^] account for the SU{3)p breaking. It was 



shown that as the symmetry breaking increases, the P always retains a larger 
binding potential than the H, and that the binding can be several tens of 
MeV. The total binding energy includes the internal kinetic energy. Because 
the c quark is massive, the kinetic energy in the P is smaller than in the H 
by about 15 MeV. This improves the prospects of the P to be bound. 

More recently, Takeuchi, Nussinov and Kubodera studied the effects 
on the Pentaquark and Hexaquark systems of instanton induced repulsive in- 
teractions for three quarks in flavor antisymmetric states. They claim in this 
framework that both Pentaquark and Hexaquark are not likely to be bound. 
Also, Zouzou and Richard reconsidered previous bag model calculations 
for the tetraquark and pentaquark. Their new calculation has weaker chro- 
momagnetic attractions at short distances and a larger bag radius for multi- 
quark states compared to ordinary hadrons. They flnd that the Pentaquark 
is unbound by 80 MeV, while the F tetraquark is unbound by 230 MeV. 
Similar conclusions for the P and H were given by Fleck et al. |T^. Riska 
and Scoccola ||18[ recently described the Pentaquark in a soliton model, using 
different parameter sets. One set gives a bound state, while another gives a 
near threshold resonance. Chow discusses the qqqcs Pentaquark in the 
framework of the binding of a heavy meson to a chiral soliton. Riska and 
Scoccola [|1^ and Oh et al. discuss the properties of a heavy qqqqc Pen- 
taquark without strangeness. Shmatikov discusses bound pentaquarks 
with a molecular type baryon- meson structure, including ND* and NB*. 

A very weakly bound -D^p deuteron-size bound state just below thresh- 
old with a structure very different from that of the strongly bound proton 
size Pentaquark might still be consistent with these recent calculations, con- 
sidering all the model uncertainties. The -D^p system does not have Pauli 
blocking and repulsive quark exchange interactions which arise in all hadron- 
hadron systems where quarks of the same flavor appear in both hadrons. 
Thus, even a comparatively weak short range interaction could produce a 
relatively large size bound state analogous to the deuteron, with a long D~y> 
tail in its wave function and a good coupling to the -D^p system. Because 
in the Pentaquark, unlike the deuteron, there is no short range repulsion, its 
structure at short distances will be quite different from that of the deuteron. 
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This component too has it's influence on the production mechanism, as dis- 
cussed in subsection 4.2. The deuteron-hke state will be stable against strong 
and electromagnetic decays. Since the D^p pair has some 50-75 MeV lower 
mass than other meson-baryon cluster components in the Pentaquark, it will 
be the dominant component in a weakly bound deuteron-like Pentaquark. 
Considering all the uncertainties in knowing the Pentaquark binding energy, 
our experimental approach is to search for both strongly and weakly bound 
Pentaquarks, as well as unbound Pentaquark resonances. 



3. Pentaquark Structure and Decay Modes 

There are different possibilities for the internal structure of observable (not 
very broad) exotic hadrons. They can be bound states or near threshold 
resonance structures of known color singlet sub-systems (e.g. AA for the H 
H or D-p for the P°). But they can have more complicated internal color 
structure, such as baryons with color octet and sextet bonds [{qqq)sc x {q(i)sc] 
and [{qqq)6c ^ {qq)6c] (see ref. [|^). We designate all such structures as direct 
five quark configurations. If color substructures are separated in space by 
centrifugal barriers, then exotic hadron resonances can have not very broad 
or even anomalously narrow decay widths, because of complicated quark 
rearrangements in the decay processes. If these exotic hadrons are bound 
strongly, they can be quasistable, with only weak decays. 

The wave function of the Pentaquark may contain two-particle clus- 
ter components, each corresponding to a pair of known color singlet par- 
ticles; and also a direct five quark component. The Pentaquark produc- 
tion mechanism and its decay modes depend on these components. The P° 
can be formed for example by the coalescence of pD~,AD ,pD*^, S+D^ + 
S°Z} ,AD ,I1'^D*~ + S°D ; or by a one-step hadronization process. Let 
us consider three color-sing let components of the P° : D'p (2907 MeV), 
L>-S+ (3058 MeV) and I^A (2981 MeV). The relative strengths of these 
components depend strongly on the binding energy, as discussed above for 
the deuteron-like Pentaquark. Pentaquark searches in progress in E791 
|2^ are based on the charged particle decay components of different 
Pentaquark decay modes: D~p (j)TT~p (B=3.5%), D~p K*^K~p 
(B_=3.3%), D-A K+TT^n'A (B=8%), WA R-n+A (B=4%) and 
D'^A K^Ti^Ti^Ti^A (B=8%). The indicated branching ratios are those 
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of the on-shell D-meson. Such Pentaquark branching ratios are plausible in 
a model where the D meson decays weakly, while the proton and A act as 
spectators. Weak decays of virtual color singlet substructures in bound states 
are possible, KD^ or for example, if their masses are smaller than the 

D~]) threshold. In other cases, there would be strong decays through quark 
rearrangement {Y,'^ D~)i,ound D~ + p, and so on. Even if the masses are 
smaller, the phase space favors decay to the lightest system. The phase space 
factor would cause the partial width for any decay mode to be smaller than 
for the on-shell decay, making the total lifetime longer. 

The decay through the direct five-quark component may open many ad- 
ditional channels, which may shorten the lifetime of the Pentaquark and 
reduce the experimental possibilities to observe it. The direct five-quark P° 
component may for example decay weakly via Cabibbo allowed or suppressed 
direct or exchange diagrams. As a result, the P° may decay into 7r~p, K~p, 
E~i^+, S+vr", S~7r+, etc. Each of these decays may have one or 

more tt+tt^ pairs, in addition to the particles shown. The observation of 
such decay channels would give important information on the Pentaquark 
internal structure. These decays may be observed, if the branching ratios are 
not too small. Depending on the decay mechanisms, the Pentaquark lifetime 
may then be shorter or longer than the Ds 467 fs lifetime. Experimental 
searches must therefore cover a range of possible lifetimes. 

We consider two possible scenarios for the decay of Pentaquark baryons. 
In the first, Pentaquarks are quasistable hadrons which can decay only 
via weak interaction, with lifetimes of the same order as other charmed 
hadrons. Possible decay modes of Pentaquarks in this scenario are 
P° — > 0pvr~, AK^TT'] P~ —>■ p(f)7T^TT~, AK'^TT~7c~ , Ti~K~^7T~] as well 
as those described in the previous paragraph. Such decays can be directly 
observed in precision vertex detectors, which are now standard devices in 
experiments with charmed or beauty particles. The use of vertex detectors 
greatly reduces the combinatorial background in the corresponding effective 
mass spectra, and makes it possible to search for exotic P baryon production 
in a wide range of the Xf observable. Searches for the tetraquark Fs (csud) 
could also look for weak decays P, such as Fg K~ K^n^n^ . 

In the second scenario, the Pentaquark baryons have large enough masses 
and are resonant states, which strongly decay with emission of secondary 
charmed particles (for example, P° pD~). The search for such strongly 
decaying exotic hadrons must also use vertex detectors for detection of 



5 



the D~. For the identification of Pentaquarks, only the study of effective 
mass spectra of secondary decay products can be used, for which there is 
a large combinatorial background. It is well known |]T], that the combi- 
natorial background is significantly reduced in the fragmentation region (at 
Xf > 0.5 — 0.6), and such kinematics is therefore strongly desirable for the 
resonance Pentaquark searches. This experimental task is more challenging 
and is crucially dependent on the Pentaquark decay width. Only narrow 
states with F < 100 MeV have a good chance to be separated from the 
background. Narrow states may arise for a variety of reasons not necessarily 
related to exotic properties, as when the phase space for decay is small. For 
example, the A(1405) is 80 MeV above the Svr threshold and has a width of 
50 MeV. The D*(2010) and the A* (2625) are about 40 MeV above their re- 
spective thresholds and their widths are less than 2 MeV. Another important 
possibility is a narrow width that may arise from the complicated internal 
color structure of an exotic hadron, and by a quark rearrangment mechanism 
in the decay processes for multiquark exotic object, leading to a colorless fi- 
nal state Another possible cause for a narrow state may be the 
reduced effective phase space as a result of OZI suppression of some decays 
for an exotic state with hidden strangeness or charm , as will be discussed 
in Section 5. 



4. Experimental Pentaquark Search 

An experimental program to search for the Pentaquark should include: 

(1) Reactions likely to produce the Pentaquark, complemented by an estimate 
of the production cross section. 

(2) Experimental signatures that allow identification of the Pentaquark. 

(3) Experiments in which the backgrounds are minimized. 

These points will be further discussed in the following subsections. 



4.1 Experimental Considerations 

All charm experiments require vertex detectors consisting of many planes of 



silicon micro-strips with thousands of channels. Fermilab E791 p3, 28 



used 23 such planes. Some of the planes are upstream of the target for beam 
tracking. These detectors allow a high efficiency and high resolution for 
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reconstruction of both primary (production) vertex and secondary (decay) 
vertex. The position resolution of the vertex detectors is typically better than 
300 microns in the beam direction. By measuring the yield of a particle as a 
function of the separation between the two vertices, the lifetime of the particle 
is obtained. This is possible as long as the hfetime is not so short, such that 
the separation of vertices becomes ambiguous. Other major components in 
charm experiments are several magnetic spectrometers with track detectors 
for track reconstruction and for momentum analysis, Cherenkov counters for 
particle identification, and electromagnetic and hadronic calorimeters. Muon 
detectors and TRD detectors for electron separation are included for studies 
of leptonic decays. The invariant mass resolution for typical charm masses 
in such spectrometers is about 10 MeV. Different spectrometers are sensitive 
to different regions of Feynman Xf values. 

In hadronic production, the charm states produced are preponderantly 
charm mesons at low Xf. The triggers for such experiments vary. In E791, 
the requirement was to ensure an interaction in the target (using signals 
from various scintillators) and a transverse energy (Et) larger than some 
threshold [|^, ^ The rest of the charm selection was done off-line. 

Future experiments are planned to obtain higher yields of charmed hadrons. 
Increased charm sensitivity can be achieved as in E781 by using higher 
integrated beam intensities, and higher efficiency detector systems. E781 
also will use a trigger condition that identifies a secondary vertex, and also 
requires positive particles with momentum greater than 15 GeV/c. This 
should enhance the high-Xj acceptance {Xf > 0.1), and give higher quality 
events. A good charm trigger can produce an enriched sample of such 
high Xf charm baryons with improved reconstruction probability because of 
kinematic focusing and lessened multiple scattering with a significantly lower 
number of events written to tape or disk. CHARM2000 experiments will also 
require charm enhancement triggers [0. The present E791 and future 
E781 and CHARM2000 experiments 0, g^, complement each other 
in their emphasis on different Xf regions, incident particle types, statistics 
and time schedules. 

High quality particle identification (PID) for the largest possible energy 
range of the outgoing particles is important for reducing backgrounds associ- 
ated with incorrect identification of tracks. In E791, two threshold Cerenkov 
detectors were used for this purpose. In E781, this will be available via ring 
imaging Cerenkov (RICH) and transition radiation detector (TRD) PID sys- 
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terns. These and other experimental techniques to reduce backgrounds are 
described in more detail in ||T^ . 



4.2 Pentaquark Production Mechanisms 

We consider possible mechanisms for P formation. For the central hadron- 
nucleus charm production at several hundred GeV/c, the elementary pro- 
cess is often associated with qq cc or gg cc transitions. The produced 
charmed quarks propagate and form mini-jets as they lose energy. Hadroniza- 
tion associated with each jet proceeds inside the nucleus, and to some extent 
also outside the nucleus; depending on the transverse momentum of the jet. 
The propagating charmed quarks may lose energy via gluon bremsstrahlung 
or through color tube formation in a string model, or by other mechanisms, 
as discussed in ref. [0] and references therein. One may form a meson, 
baryon, or Pentaquark, according to the probability for the charmed quarks 
to join together with appropriate quarks and antiquarks in the developing 
color field. One can estimate Pentaquark production cross sections via one- 
step and also two-step hadronization. All such estimates are very rough. 
Our aim is to account for major ingredients in estimating the cross section, 
and to give a conservative range of values. For one-step hadronization, the 
c joins directly to the other quarks. The one-step is the usual mechanism 
for meson and baryon formation. For two-step, the first involves meson and 
baryon hadronization, while the second involves meson-baryon coalescence. 

We first consider estimates for the central production cross section as- 
suming a meson-baryon coalescence mechanism, expected to be the main 
mechanism for production through the long-range (deuteron-like) compo- 
nent of the Pentaquark wave function. We make a crude estimate relative to 
the D~, an anticharmed-strange meson (cs). The weakly bound P (deuteron 
type structure) can be produced for example by coalescence of a proton or a 
neutron with a D~, analogous to the production of a deuteron by coalescence 
of a neutron and a proton. The data give roughly 10~^ for the a{d)/a{p) 
production ratio. This ratio can also be applied to a{P) / a[D~) production. 
The reason is that in both cases, the same mass (nucleon mass) is added 
to the reference particle (proton or D7), in order to form a weakly bound 
deuteron-like state. This ratio is very sensitive to the pentaquark binding 
energy, and may be substantally different for a larger value. 

We now consider the one-step hadronization of a Pentaquark, expected to 
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be the main mechanism for the production through the short-range compo- 
nent of the Pentaquark wave function. We rely here on an empirical formula 
which reasonably describes the production cross section of a mass M hadron 
in central collisions. The transverse momentum distribution at not too large 
Pi follows a form given as ||36| : 



da /dpi ~ exp{-C^jM^ + (4) 
where C is roughly a universal constant ~ 5 - 6 (GeV)~^. The exponential 



(Boltzmann) dependence on the transverse energy Et = \J + pi has in- 
spired speculation that particle production is thermal, at a temperature 
~ 160 MeV We assume that this equation is applicable to Pentaquark 
production. To illustrate the universality of C, we evaluate it for a few cases. 
For Ac and empirical fits to data give exp(-6j9j), with b=l.l GeV^^and 
b=2.0 GeV-2, respectively |^ Hf. With C ^ 2 ■ 6 ■ M, this corresponds 



to C ~ 5.0 GeV ^ for and C ~ 5.3 GeV ^ for For inclusive pion 
production, experiment gives exp(-6pi) with b = 6 GeV^^ p9i; and C ~ b, 



since the pion mass is small. Therefore, C= 5-6 GeV^^ is valid for Ac, S° 
hyperon, and pion production. After integrating over p^, we estimate the 
ratio as: 

a{P)/a{D-) ~ exp[-h[M{P) - M(D;)]] ~ lO'^. (5) 

In applying Eq. 5 to Pentaquark production, we assume that the suppression 
of cross section for the heavy P as compared to the light D~ is due only to 
the increased mass of P. The formula ignores dynamical considerations, such 
as the particular one-step or two-step hadronization processes, or the size of 
the P. It also does not account for threshold effects. Consequently, there may 
be large uncertainties in its application. 

Both reaction mechanisms described above can contribute to the pro- 
duction cross section, which is estimated in the range of a{P) / a{D~) ~ 
10~^ — 10~^. In actual measurements, the product a ■ B for a particular 
decay mode is measured, where B is the branching ratio for that mode. If 
a Pentaquark peak is not observed, assumptions on the values of B and on 
the P lifetime may be necessary in order to set limits to the Pentaquark 
production cross section. 

Another approach for the cr(P) cross section evaluation is based on a 
comparison of Pentaquark and charmed-strange baryon H° (csd) production 
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reactions. This method was previously apphed in p4| , p5| , p6| . For example, 
with a high energy hyperon beam, one may compare P° production (uudcs) 
to (esc?) baryon production. Processes with charm baryon production 
involve cc pair creation, followed by c or c hadronization in the final states. 
There should therefore not be much difference for c or c fusion, so that one 
may compare P° to H° production. The production of P° is also different 
from S° due to the fusion of an additional uu quark pair. This would cause 
a reduction factor R for the cross sections: 

as(P°) = > 2.5 x 10-^^(2°). (6) 

This factor R > 2.5 x 10^^ can be estimated from the data on other processes 
with additional quark fusion in the particles under study. For example, the 
relative yields d/p ~ 10~^, •^He/p ~ 10~^ ^4, give information 



on the fusion probability of 3 and 6 additional quarks. It is then possible to 
obtain an average reduction factor ~ 5 x 10~^ per fusion for each additional 
quark, yielding the value R > 2.5 x 10^^ in Eq. 6. The inequality arises 
because antinuclei are loosely bound systems whose yields implicitly include 
this factor. 



4.3 Pentaquark Decay Signatures 

(1) Mass, Width and Decay Modes 

Searches for the Pentaquark are easiest via modes having all final decay 
particles charged. With all charged particles detected, the invariant mass of 
the system can be determined with high resolution. One signature of the 
Pentaquark is a peak in the invariant mass spectrum somewhat lower than 
2907 MeV if the system is bound; and above if it is a resonance. The position 
of the peak should be the same for several decay modes. It's width should 
be determined by the experimental resolution if it is bound, and broader if 
it is a resonance. 

The selection of the decay modes to be studied is made primarily by 
considering detection efficiency and expected branching ratios. Since the 
D~p system is the lightest, it is expected to be preferred from phase space 
arguments. Also, two of it's decay modes have four charged particles in 
the final state (e.g. (pir-p, ^ K+R-; K*^K-p, K*^ K+ir-). This 



signature was implemented in E791 p4|, BH]. First, two distinct vertices were 
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identified, a production vertex and a decay vertex. From the decay vertex, 
four tracks were identified and associated with K^K^n~p. By reconstructing 
the invariant mass of the K^K~ pair, only mass events were accepted 
and the invariant mass of all four particles was reconstructed. A peak in 
the resulting spectrum will be one of the identifying characteristics of the 
Pentaquark. 

(2) One General Signature - A Spectator Baryon: 

We first note a striking signature for Pentaquark decay which may be 
useful for discrimination against background. This signature is predicted 
by both of two very different Pentaquark models: (1) a loosely-bound D~p 
deuteron-like state and (2) a strongly-bound five-quark state. Both models 
predict decay modes into a baryon and two or more mesons, in which the 
three quarks in the baryon are spectators in the decay process and remain in 
the final state with a low momentum which is just the fermi momentum of 
the initial bound state. 

That the baryon is a spectator is obvious in the deuteron model, in which 
the decay is described as an off-shell D~ decaying with a nucleon spectator. 
In the five-quark model, a similar situation arises in the commonly used 
spectator model with factorization. Here, the charmed antiquark decays 
into a strange antiquark by emission of a W~ which then creates a quark- 
antiquark, which hadronizes into mesons. The strange antiquark combines 
with one of the four spectator quarks to form one or more mesons, while the 
three remaining spectator quarks combine into a baryon. 

In both cases, it seems that the final state should show a low-momentum 
baryon in the center-of-mass system of the Pentaquark and the invariant mass 
spectrum of the remaining mesons peaked at the high end near the kinematic 
limit. Thus in the particular cases of the p4>n~, K*°K~p and AK^n^ decay 
modes, the (j)7T~, K*°K~ and X+tt" invariant mass distributions respectively 
should show this peaking near the kinematic limit. 

Note that in the particular case of the pcpn' decay mode, a low momentum 
proton in the center of mass system means that the vr^ and arc back to 
back with the same momentum and therefore that the pion carries off most 
of the available energy. Thus one might reduce background with a cut that 
eliminates all pions with low momentum in the center of mass. 

(3) Some Model-Dependent Branching Ratio Predictions: 

The 07r~j9 decay mode is the most convenient for a search, since the (f) 
signal is so striking. We now examine the lowest order predictions from the 
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two extreme models for the branching ratios of other modes relative to (f)7r~p. 

In experiments sensitive only to charged particles the (pn^p decay mode 
is observed in the four-prong final state K^K~T^~p. The K*°K'~p decay 
mode is also observable in this same four prong final state. The K*°K~p 
decay mode arises naturally in the deuteron model, since the K*°K~ decay 
is observed for decays with a comparable branching ratio to (j)TX~ . In 
this model, the ratio of the two decays is predicted from observed decay 
branching ratios with phase space corrections. However, the K*°K~p decay 
mode does not occur in the five quark spectator model, where the spectator 
strange quark can only combine with the s produced by the charm decay 
to make a or with two spectator nonstrange quarks to make a hyperon. 
Comparing the two decays thus tests the decay model. 

The XttA and K*t^K decay modes arise naturally in the five quark specta- 
tor model, or in a moderately bound DA Pentaquark. However, they should 
not be expected in a very weakly bound deuteron model with mainly a D~y> 
structure. In that case, the decays into mesons containing one strange 
quark-antiquark pair and the baryon spectator has no strangeness. 

(4) Angular Momentum Constraints and Angular Distributions for P De- 
cays: 

We can give a model-independent prediction. The Pentquark has spin 
1/2 and this total angular momentum is conserved in the decay. Since the 
production process is a strong interaction which conserves parity, the Pen- 
taquark will not be produced with longitudinal polarization. Its polarization 
in the beam direction must also vanish. Therefore, the angular distribution 
in the center-of-mass system of the Pentaquark must be isotropic for the mo- 
mentum of any final state particle in any decay mode with respect to either 
the incident beam direction or the direction of the total momentum of the 
Pentaquark. The background does not necessarily have these constraints. 

We also give a model-dependent prediction. We first consider the 
deuteron model. The D~ has spin zero, and spin is preserved in the de- 
cay. Thus, in the center of mass frame of all the D~ decay products, the 
angle between the proton momentum and the momentum of any particle 
emitted in the Dj decay must have an isotropic angular distribution. 

A further prediction is obtainable for the case of a vector-pseudoscalar 
decay mode of the D~\ e.g. 07r~ or K*^K~ . The vector meson must be 
emitted with zero helicity in the rest frame of the D~ . The zero helicity 
can be seen in the 07r decay by measuring the angle Oktx between the kaon 
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momenta in the rest frame and the pion momentum. The prediction is 
to have a cos^ Ok-k distribution. By contrast, the five-quark model for the 
Pentaquark favors hehcity one over hehcity zero for the vector meson by just 
the 2:1 ratio needed to give an isotropic distribution in Ok-k- Here again the 
background does not necessarily have these constraints. 



4.4 Pentaquark Expected Yield 

We proceed with count rate estimations for the expected yields of Pentaquark 
baryons. Given the need to search for both quasistable and resonant Pen- 
taquark baryons in different regions of X/, we give production cross sections 
ioiXf > and in the fragmentation region > 0.5, where one expects an 
improved signal to background ratio. We use several production models for 
Pentaquarks. The different predictions of these models reflect the uncertain- 
ties in our expectations. We also assume that quasistable P° baryons would 
be reconstructed in selected visible weak (w) decay modes with a combined 
effective branching ratio B^: 

= B[P^ ^ p(f)Tr- + p(f)Tr-TT+TT- + pK'^*K-] 
5^ = ~ 0.05 ; based on (p ^ K+R-and ^ K+n~. (7) 

In practice, one looks for the visible (charged-particle) decay modes, such 
as: [D-p] <pTi-p K+K-n-p (B=1.8%), [D-p] (pTi-ji-Ti+p 
K+K-n-n-n+p (B=0.9%), and [D-p] K*^K-p n-R+R-p 

(B=2.1%). The total branching ratio of these last three decay modes is 
4.8%. Here we assume a [D~p] bound state model for the quasistable P 
pentaquark. Our estimate By^ =~ 0.05 is somewhat lower than the total 
visible (charged particle) decay modes of a virtual D~ for P° {pD^)mrtuai 
dissociation, where B{Dj) (visible) ~ 0.08 The 0.08 value includes 



possible non-resonant charged particle decays, such as D~ — > pR^R~n~ or 
D~ —>■ p'K~^TT~7r~. Our objective is a search that will have improved signal to 
background, by concentrating on final states having resonances, such as the 
narrow or the K*(890). 

For strongly-decaying P baryon resonances, we assume the same produc- 
tion cross sections as for a quasi-stable Pentaquark. A reasonably small 
decay width can be obtained in principle in the model with direct five-quark 
configurations (color octet or sextet bonds, as described in Section 3 and 
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Ref. |2^). For detection of these baryons, one may hunt for the decay mode 
P° —* pD^ , assigned here a 100% branching ratio for this resonance possi- 
bihty. This value is clearly an upper limit, and may be reduced by other 
strong decay channels. We for example do not take into account other pos- 
sible totally reconstructable decay modes, such as P° pD^ {visible)TT^TT^ , 
— > D~ {visible) . Thus, we estimate for strongly decaying pen- 
taquarks, the branching ratio: 

B,[P°]{visible) = = B[P° pD;] ■ B[D; visible] ~ 0.05. (8) 

One may also search for P Pentaquark strong decays P D^A; D^I1~ 
(resonances). 

Both weak and strong decay modes coming from the D~p and the D^A 
components of the P are currently being studied in E791, where the data 
were taken with a 500 GeV vr" beam. Analysis of a part of the E791 data 
already yielded a preliminary upper limit at 90% confidence level: 

for quasistable Pentaquark production (not including systematic uncertain- 
ties) [Q. In the analysis, it was assumed that the Pentaquark has the same 
lifetime as a real D~ and a mass of 2.75 GeV. It was assumed that its produc- 
tion characteristics are the same as other charmed baryons. This limit was 
based on a part of the data and measured D~ yield. With the full data sam- 
ple and more decay modes analyzed, several Pentaquarks may be observed 
if the cross section is in the range estimated in the previous section, or else 
the limit may be further lowered. 

For the planned E781 and CHARM2000, when both use baryon beams, 
we rely on previous measurements done with similar beams. The S~ hyperon 
beam should be good tool for the search for strange-charmed(anticharmed) 
open exotic hadrons, such as the Pentaquark baryons P° and P~, or 
tetraquark meson P,. Hyperon beams are the purest high energy beams 
containing strange valence quarks. Thus, the inclusive reactions: 

S-(rfrfs) + ^ P^{cuuds) + X 
P~ {cudds) + X 

F^{csud)+X (10) 
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should be favorable for the production of strange-charmed hadrons at large 
Xj, since they benefit from strange quark sharing between primary and sec- 
ondary particles. 

We use different methods to estimate the expected cross sections for Pen- 
taquark production in I1~N interactions for the hyperon beam of Fermilab 
with momentum 650 GeV/c (E781). First, following section 4.2, we take 
(Js-(-P°)/o-s-(-D7) ^ 10~^ — 10~^. The aY;-{D~) cross section is estimated 
two ways. A neutron beam measurement with En — 600 GeV gave 
aniD;) ■ B{D- <pTi-) ~ 0.76/i6/A^ for 0.05 < X/ < 0.35. 

From the quark structure of and the neutron, one may expect that 
o'sl-D^) is greater than an{D~). Assuming conservatively that they are 
equal, and assuming also the Xf dependence da{D~) / dXf ~ (1 — X/)^ 
for the baryon primaries (see ||3^), one obtains crY.{D~)\xf>o — 40 jjb/N. 



This is an exceptionally large value, considering that the entire hadronic 
(cq) production cross section is expected |Q to be near 20-30 /ib/N. A 
more conservative estimate of this cross section is based on the assumption 
aY,{D~) ~ ap{D~) ~ ap{D~^) from the quark structure of the projectiles 
and produced mesons. The data of the EGS experiment for Xf > with a 
Pp = 400 GeV/c momentum proton beam give (Tp(i5^) ~ cXp{D^) ~ 3 ^b/N 



45|] . Based on the energy dependence of charm production data , extrap- 
olation to 650 GeV/c yields for E781 the estimate cr-£{D~)\xf>o — 5 jjb/N. 
This smaller value is used for the estimate of P° exotic baryon yields 
in E781 (see Table 1). With a{P'^)/a{D~) ~ lO^^ - IQ-^, one obtains 
(t(P°) ^50-5 nh/N. 

There are no direct data for the production with baryon beams, and 
only very poor data for (csu) production [|^. The recent WA89 experi- 
ment used a S~ hyperon beam at CERN, with momentum Ps-=330 GeV/c. 
The cross section measured was: crE-(A+)|xj.>o.2 = (9.3 ± 4.3 ± 2.5) fJ^b/N, 
with da^/dXf ~ (1 - X/)^ [0]. Extrapolating to X/ > (factor ~ 3) 
and to Ps- =650 GeV/c (factor ~ 1.5), one obtains as- (A+)|x^>o — (30± 
14 ± 8) fib/N for E781, where the uncertainties are only from scaling those 
of Ref. The value 30 /ib/N again seems unreasonably large, which 

may be due to the extrapolation procedure or to the experimental uncer- 
tainties in o"e-(A+)|x^>o.2- Thus, we use the more conservative estimate 
crs-(A+)|x^,>o ~ 5 fib/N for Ps-=650 Gev/c, close to the minimum value 
given by the error bars. From the quark composition of charmed baryons, it 
is reasonable to conclude that (Ts(S°) > crs(H+) or crs(A+) > 5 fib/N; 
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and to obtain from Eq. (6) the estimate: 



MP')\x,>o > 13nb/N. (11) 

This value gives Pentaquark yields in E781 at about 25% of the upper values 
shown in Table 1. 

We consider also the expected E781 efficiency for Pentaquark detection, 
by comparison to estimated efficiencies for cqq decays. These include 
a tracking efficiency of 96% per track, a trigger efficiency averaged over Xf 
of roughly 18%, and a signal reconstruction efficiency of roughly 50%. The 
E781 Monte Carlo simulations ||2^ gave an average global efficiency of ~ 
8%, by considering the ~ 200 fs lifetime decay A+ pK^n^, and the ~ 



350 fs lifetime decay S^tt+tt^. The charm baryons were assumed p9 

to be produced with a cross section of the form da/dXf = (1 — X/)^'^, an 
assumption which is built into the estimation of the trigger efficiency. For 
heavier Pentaquark production, it is likely that this distribution would shift 
to lower Xf (corresponding to an exponent greater than 4.2), and this would 
also reduce the efficiency. The value 8% is for reconstruction of the relatively 
strong signals from cqq charm baryon decays. The reconstruction efficiency 
may be lower for Pentaquark events. For the weaker Pentaquark signal, 
tighter analysis cuts with resulting lower efficiencies may be required in order 
to achieve the optimum signal to noise ratio. For lifetimes smaller than about 
60 fs, which is possible for the Pentaquark, the trigger efficiency would also 



be significantly reduced p9[. Considering all the unknown variables, the 
final experimental statistics may then be significantly lower than the upper 
value estimated here, using the 8% global efficiency. From Table 1, one sees 
that for quasistable Pentaquark baryons, the maximum expected statistics 
in E781 is 280-2800 events. This may be adequate for their observation, if 
they really exist. 

From da/dXf ~ (1 — XjY for charmed baryons |3^, we estimate the 
Pentaquark production cross section in the Xj > 0.5 fragmentation region 
(quasi-stable or resonant) as ad ~ 0.03 ■ cr{Xf > 0); as given in Table 1. The 
fragmentation region should be most effective for reducing the combinatorial 
background in both quasi-stable and resonance Pentaquark searches. In the 
latter case, one studies a strong decay into D~p, if the P is a short-lived 
resonance with mass M{P^) > M{D^) + M{N). For this strong decay, 
the proton and D~ come from primary vertex, and the D~ decay forms 
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the secondary vertex. The expected maximum numbers of events in the 
fragmentation region is quite hmited (from 28 to 280 events, from Table 1), 
and there are moderate chances for P*^ observation as a peak in the mass 
spectrum of M{pD~). 

It is possible that different mechanisms for charm production contribute 
in different Xf regions. For example, there is evidence for leading production 
of charmed hadrons in WA89 and FNAL E769 ||] and E791 ||. One can 



also consider a diffractive mechanism for Pentaquark production. Brodsky 
and Vogt [^, ^ suggested that there may be significant intrinsic charm (IC) 
cc components in hadron wave functions. The Hoffmann and Moore analysis 
of charm production in deep inelastic electron scattering yields 0.3% IC 
probability in the proton. A recent reanalysis of the EMC charm production 



data was carried out by Harris, Smith, and Vogt ||5^ . Their improved analysis 
found that an IC component is still needed to fit the EMC data, with a value 
indicated for the proton of (1.0 ± 0.6)%. The most probable IC state occurs 
when the constituents have the smallest invariant mass. In the rest system, 
this happens when the constituents are relatively at rest. In a boosted frame, 
this configuration corresponds to all constituents having the same velocity 



and rapidity 55]. When a IC state is freed in a hadronic collision, the 
charm quarks should have approximately the same velocity as the valence 
quarks. They can then easily coalesce into charmed hadrons Q ^ Q 
and produce leading particle correlations at large Xf. The IC component 
in an incident or proton can then lead to large Xf P° production. One 
may expect that P° production will be predominantly central for reaction 
mechanisms other than IC. Intrinsic charm Pentaquark production, with its 
expected high X/ distribution, would therefore be especially attractive. 

For E781 and CHARM2000 one can study the diffractive pair production 
reactions ([0, see also 0) S^+iV ^ {P-D^)+N and p+N {P°D+)+N 
(or even the coherent reactions of these types on nuclei), with possible D^, 

tag or without such tag. For the diffractive pair production cross section, 
one can compare to the diffractive cross section for the reaction p + N ^ 
{KK^) + iV at 70 GeV, which is more than 1 iib after subtraction of isobar 



contributions [^]. The ratio of diffractive cross sections ad was estimated 



ad{P^D+)[QmGeV]|ad{^K+)[^QGeV] ~ (m./mj^ . R ^ 2. x 10"^ (12) 
Here the reduction factor R that accounts for the fusion in the P of an extra 
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two quarks is R=2.5 x 10~^ from Eq. 6. The factor {ms/mc)'^ accounts for 
the relative probabihty to produce charmed quark versus strange quark pairs. 
From the ratio of constituent quark masses [^], (m<j/mc)^ ~ 8. x 10~^. In 
Eq. 12, we do not exphcitly show a kinematic factor, since its value |l|, ^ 



is close to unity for the reactions shown. This factor accounts for the energy 
and mass dependence of the cross sections, for diffract ive-like production 
of different final state masses with different beam energies. The increased 
Pentaquark cross section expected from energy extrapolation is offset by a 
reduction for making the heavier mass Pentaquark final state. We then find: 

ad{P^D+)[QOOGeV] ~ 1 /i6 x 2. x 10"^ ~ 0.2 nb. (13) 

We assume conservatively that -B„,=0.05 for such diffractive production 
searches, although with the smaller backgrounds expected at high Xf, one 
could possibly search for all visible decays (5^=0.08). From the projected 
Pentaquark charm sensitivity shown in Table 1 for E781 (and roughly ten 
times higher for CERN CHEOPS it may be possible to observe the 

diffractive production process of Pentaquarks at the level of several hundred 
events. Future experiments in the spirit of the CHARM2000 workshop, with 
higher charm sensitivity than E781, should have improved chances to ob- 
served Pentaquarks. 



5. Heavy Baryons with Hidden Charm 

In recent years, several candidates were reported for baryon states with un- 
usual properties (narrow decay widths, large branching ratios for the decays 
with strange particles). There are candidates for cryptoexotic baryons with 
hidden strangeness =\ qqqss>, were q = u or d quarks (see [|I|, |l^, ^ 
and references therein). Further searches for nonstrange and cryptoexotic 
baryons with hidden strangeness are planned at IHEP (70 GeV proton beam) 
with the SPHINX spectrometer and in FNAL E781 |ET|. Although the 
existence of such a baryon is not yet confirmed, this suggestion raise the 
question of the possible existence of heavy cryptoexotic baryons with hidden 
charm B^ =\ qqqcc>. These exotic Pentaquark baryons can either be direct 
five-quark states or (p — N and — N bound states (resonances). The lat- 
ter may form at low relative velocities of the meson and baryon consituents. 
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because of the strong attractive QCD Van der Waals interaction [B8[. The in 



trinsic charm which was discussed in Section 4.4 may also be relevant for the 
production of the baryon. If M{B^) <M(r/J+M(p) ^ 3.9 GeV, the 
decays would be OZI suppressed and the width of this cryptoexotic baryon 
would be quite narrow. If M{B^) >4.3 GeV, there would be OZI allowed 
decays B^ —>■ p+ J/ip; + Z)°, etc. Because of a complicated internal color 
structure of this baryon (see Introduction), one can expect a narrow decay 
width (< 100 MeV). 

To search for such B^ states, it was proposed to use the diffrac- 
tive production reaction p + N B^ + N with possible decays of B^ 
baryons B^ ^ p + {J/ip)virt_ ^ p + or B^ ^ p + {'nc)virt 

p+ {K~^K^'K^'n-^;2'n-^27r~; KKn;r]7nr). For these processes, aBe was es- 
timated as 0.2 - 0.5 nb/N [^]. For an appropriately designed experiment, 
with a charm sensitivity close to 10^ events/(nb/N) of effective cross section, 
it should be possible to observe several thousand such events. 



6. Conclusions 

We described the expected properties of Pentaquarks. Possibilities for en- 
hancing the signal over background in Pentaquark searches were investigated. 
General model-independent predictions were presented as well as those from 
two models: a loosely bound D~N "deuteron" and a strongly-bound five- 
quark model. While the current E791 may have marginal sensitivity, future 
experiments with more than 10^ reconstructed charmed baryon events should 
have enough sensitivity to determine whether or not the Pentaquark exists. 
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Table 1: Projected Pentaquark Yields 



Model for estimation 


a{pO)/a{DJ) ~ 10-2 - 10-3 


of 




~ 5 ^lb/N 


at 


- = 650 GeV 


FNAL E781 


pO 


'^E-(-P")lx,>o 


~ 50 - 5 nb/N 




0.05 


quasi- 
stable; 


£1 


< 0.08 


Xf>0 




< 0.2 - 0.02 nb/N 






< 2800 - 280 events 


pO 




~ 1.5 - 0.15 nb/N 


Bs or Bu, 


0.05 




£2 


< 0.25 


Xf > 0.5 




< 0.02 - 0.002 nb/N 




Ns or 


< 280-28 events 



Notes to Table 1: 



1. Here u-^- C^'") = ""{^ -|- Af — > P" + X), etc. Bw is the effective branching ratio for visible weak 
decays of quasistable exotic Pentaquark baryons: 

Bi„ = B[P° pcpTT- +p4>w--rT+TT- +pK°*K-] ~ 0.05, as described in the text. 

Bs is the effective branching ratio for visible decays of a Pentaquark exotic baryon resonance: 

Bs = B[P° pDJ] ■ B[DJ visible] ~ (1.0) ■ (0.05) ~ 0.05. 

o-(pO)^jf = o-j,_ (pO)._B £, where B is either B„ or Bs. We use £i < 0.08 for the average efficiency 
for the reconstructed charm events i n F, 781 for X f > 0; ei < (trigger efficiency)- (reconstruction 
efficiency) < (0.18)- (0.45) = 0.08 psj]. The efficiency depends on the lifetime, and also on the 
particular decay mode observed, and also on the average value for Pentaquark events. We use 
£2 < (0.55) ■ (0.45) = 0.25 for the same efficiency in the fragmentation region {Xf > 0.5); since 
the trigger efficiency is higher for higher Xf particles. We estimate the Pentaquark production 
cross section in the fragmentation region (quasi-stable or resonant) as ~ 0.03 ■ (^{Xf > 0). N^j 
is the number of weak decay events for quasistable P". Ns is the number of decay events for a 
resonance P^ . Th e estimated events are extrapolations for a planned Fermilab E781 1996-97 50 
week data run tesl]: 1.3 X 10^^ interactions in the target, an estimated 2.8x10* charm events, and 
a sensitivity of roughly 1.4 X 10'' events/ (nb/N) of effective charm cross section. The values cited 
for N™ and Ns are only projected upper limits, as described in the text. 
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